Variants at the ASIP locus contribute to coat color darkening in Nellore cattle.

BACKGROUND: Nellore cattle (Bos indicus) are well-known for their adaptation to warm and humid environments. Hair length and coat color may impact heat tolerance. The Nellore breed has been strongly selected for white coat, but bulls generally exhibit darker hair ranging from light grey to black on the head, neck, hump, and knees. Given the potential contribution of coat color variation to the adaptation of cattle populations to tropical and sub-tropical environments, our aim was to map positional and functional candidate genetic variants associated with darkness of hair coat (DHC) in Nellore bulls. RESULTS: We performed a genome-wide association study (GWAS) for DHC using data from 432 Nellore bulls that were genotyped for more than 777 k single nucleotide polymorphism (SNP) markers. A single major association signal was detected in the vicinity of the agouti signaling protein gene (ASIP). The analysis of whole-genome sequence (WGS) data from 21 bulls revealed functional variants that are associated with DHC, including a structural rearrangement involving ASIP (ASIP-SV1). We further characterized this structural variant using Oxford Nanopore sequencing data from 13 Australian Brahman heifers, which share ancestry with Nellore cattle; we found that this variant originates from a 1155-bp deletion followed by an insertion of a transposable element of more than 150 bp that may impact the recruitment of ASIP non-coding exons. CONCLUSIONS: Our results indicate that the variant ASIP sequence causes darker coat pigmentation on specific parts of the body, most likely through a decreased expression of ASIP and consequently an increased production of eumelanin.

Chromosomal segments may explain the antibody response cooperation for canine leishmaniasis pathogenesis.

Visceral leishmaniasis (VL) is marked by hyperactivation of a humoral response secreting high quantity of immunoglobulins (Igs) that are inaccessible to intracellular parasites. Here we investigated the contributions of the antibody response to the canine leishmaniasis pathogenesis. Using correlation and genome-wide association analysis, we investigated the relationship of anti-Leishmania infantum immunoglobulin classes levels with parasite burden, clinical response, renal/hepatic biochemical, and oxidative stress markers in dogs from endemic areas of VL. Immunoglobulin G (IgG) and IgA were positively correlated with parasite burden on lymph node and blood. Increased IgG, IgA and IgE levels were associated with severe canine leishmaniasis (CanL) whereas IgM was elevated in uninfected exposed dogs. Correlations of IgM, IgG and IgA with creatinine, urea, AST and ALT levels in the serum were suggested an involvement of those Igs with renal and hepatic changes. The correlogram of oxidative radicals and antioxidants revealed a likely relationship of IgM, IgG and IgA with oxidative stress and lipid peroxidation in the blood, suggested as mechanisms mediating tissue damage and CanL worsening. The gene mapping on chromosomal segments associated with the quantitative variation of immunoglobulin classes identified genetic signatures involved with reactive oxygen species generation, phagolysosome maturation and rupture, free iron availability, Th1/Th2 differenciation and, immunoglobulin clearance. The findings demonstrated the roles of the antibody response as resistance or susceptibility markers and mediators of CanL pathogenesis. In addition we pinpointed candidate genes as potential targets for the therapy against the damage caused by exacerbated antibody response and parasitism in VL.

Misidentification of runs of homozygosity islands in cattle caused by interference with copy number variation or large intermarker distances.

BACKGROUND: Runs of homozygosity (ROH) islands are stretches of homozygous sequence in the genome of a large proportion of individuals in a population. Algorithms for the detection of ROH depend on the similarity of haplotypes. Coverage gaps and copy number variants (CNV) may result in incorrect identification of such similarity, leading to the detection of ROH islands where none exists. Misidentified hemizygous regions will also appear as homozygous based on sequence variation alone. Our aim was to identify ROH islands influenced by marker coverage gaps or CNV, using Illumina BovineHD BeadChip (777 K) single nucleotide polymorphism (SNP) data for Austrian Brown Swiss, Tyrol Grey and Pinzgauer cattle. METHODS: ROH were detected using clustering, and ROH islands were determined from population inbreeding levels for each marker. CNV were detected using a multivariate copy number analysis method and a hidden Markov model. SNP coverage gaps were defined as genomic regions with intermarker distances on average longer than 9.24 kb. ROH islands that overlapped CNV regions (CNVR) or SNP coverage gaps were considered as potential artefacts. Permutation tests were used to determine if overlaps between CNVR with copy losses and ROH islands were due to chance. Diversity of the haplotypes in the ROH islands was assessed by haplotype analyses. RESULTS: In Brown Swiss, Tyrol Grey and Pinzgauer, we identified 13, 22, and 24 ROH islands covering 26.6, 389.0 and 35.8 Mb, respectively, and we detected 30, 50 and 71 CNVR derived from CNV by using both algorithms, respectively. Overlaps between ROH islands, CNVR or coverage gaps occurred for 7, 14 and 16 ROH islands, respectively. About 37, 44 and 52% of the ROH islands coverage in Brown Swiss, Tyrol Grey and Pinzgauer, respectively, were affected by copy loss. Intersections between ROH islands and CNVR were small, but significantly larger compared to ROH islands at random locations across the genome, implying an association between ROH islands and CNVR. Haplotype diversity for reliable ROH islands was lower than for ROH islands that intersected with copy loss CNVR. CONCLUSIONS: Our findings show that a significant proportion of the ROH islands in the bovine genome are artefacts due to CNV or SNP coverage gaps.

Canine leishmaniasis: Genome-wide analysis and antibody response to Lutzomyia longipalpis saliva.

The anti-inflammatory properties of sand fly saliva favor the establishment of the Leishmania infantum infection. In contrast, an antibody response against Lutzomyia longipalpis saliva is often associated with a protective cell-mediated response against canine visceral leishmaniasis. Genetic studies may demonstrate to what extent the ability to secrete anti-saliva antibodies depends on genetic or environmental factors. However, the genetic basis of canine antibody response against sand fly saliva has not been assessed. The aim of this study was to identify chromosomal regions associated with the anti-Lu. longipalpis salivary IgG response in 189 dogs resident in endemic areas in order to provide information for prophylactic strategies. Dogs were classified into five groups based on serological and parasitological diagnosis and clinical evaluation. Anti-salivary gland homogenate (SGH) IgG levels were assessed by Enzyme-Linked Immunosorbent Assay (ELISA). Genomic DNA was isolated from blood samples and genotyped using a SNP chip with 173,662 single nucleotide polymorphism (SNP) markers. The following linear regression model was fitted: IgG level = mean + origin + sex + age + use of a repellent collar, and the residuals were assumed as pseudo-phenotypes for the association test between phenotypes and genotypes (GWA). A component of variance model that takes into account polygenic and sample structure effects (EMMAX) was employed for GWA. Phenotypic findings indicated that anti-SGH IgG levels remained higher in exposed and subclinically infected dogs than in severely diseased dogs even in regression model residuals. Five associated markers were identified on chromosomes 2, 20 and 31. The mapped genes included CD180 (RP105) and MITF related to the rapid activation of B lymphocytes and differentiation into antibody-secreting plasma cells. The findings pointed to chromosomal segments useful for functional confirmation studies and a search for adjuvant molecules of the anti-saliva response.

Pleiotropic Genes Affecting Carcass Traits in Bos indicus (Nellore) Cattle Are Modulators of Growth.

Two complementary methods, namely Multi-Trait Meta-Analysis and Versatile Gene-Based Test for Genome-wide Association Studies (VEGAS), were used to identify putative pleiotropic genes affecting carcass traits in Bos indicus (Nellore) cattle. The genotypic data comprised over 777,000 single-nucleotide polymorphism markers scored in 995 bulls, and the phenotypic data included deregressed breeding values (dEBV) for weight measurements at birth, weaning and yearling, as well visual scores taken at weaning and yearling for carcass finishing precocity, conformation and muscling. Both analyses pointed to the pleomorphic adenoma gene 1 (PLAG1) as a major pleiotropic gene. VEGAS analysis revealed 224 additional candidates. From these, 57 participated, together with PLAG1, in a network involved in the modulation of the function and expression of IGF1 (insulin like growth factor 1), IGF2 (insulin like growth factor 2), GH1 (growth hormone 1), IGF1R (insulin like growth factor 1 receptor) and GHR (growth hormone receptor), suggesting that those pleiotropic genes operate as satellite regulators of the growth pathway.